Electron spin resonance shifts in S = 1 antiferromagnetic chains 
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We discuss electron spin resonance (ESR) shifts in spin-1 Heisenberg antiferromagnetic chains with 
a weak single- ion anisotropy based on several effective field theories, the 0(3) nonlinear sigma model 
(NLSM) in the Haldane phase, free fermion theories around the lower and the upper critical fields. 
In the 0(3) NLSM, the single- ion anisotropy corresponds to a composite operator which creates two 
magnons at the same time and position. Therefore, even inside a parameter range where free magnon 
approximation is valid, we have to take interactions among magnons into account. Though the 0(3) 
NLSM is only valid in the Haldane phase, an appropriate translation of Faddeev-Zamolodchikov 
operators of the 0(3) NLSM to fermion operators enables one to treat ESR shifts near the lower 
critical field in a similar manner to discussions in Haldane phase. We present that our theory 
gives quantitative agreements with recent ESR experimental results on an spin-1 chain compounds 
NDMAP. 

PACS numbers: 76.20.+q, 75.10.Jm, 75.30.Gw 
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I. INTRODUCTION 

Quantum phase transition has been studied for a long 
time. In quantum magnetism, the magnetic field is the 
most familiar parameter to cause quantum phase transi- 
tions. An 5 = 1 Heisenberg antiferromagnetic (HAF) 
chain and an S = 1/2 two-leg HAF ladder are typi- 
cal examples of one-dimensional quantum spin systems 
which show the quantum phase transition by the mag- 
netic field. These systems have the unique ground state 
separated from excited states by a finite excitation gap at 
zero field. As the magnetic field is gradually applied, the 
excitation gap is going to vanish*^ After the collapse of 
the excitation gap, the system enters into a field-induced 
critical phase. The field-induced critical phase lies in a 
range H c \ < H < H C 2- Here H c i and H C 2 are called as 
a lower and an upper critical field. For H < H c i, the 
system is in the gapped phase. And for H C 2 < H, the 
system is in another gapped phase where the spins are 
fully polarized. H C 2 is also called as a saturation field. 

The quantum phase transitions at H = H c \ and 
H = H C 2 bring about reconstructions of the excitation 
spectrum. Especially, dynamical properties of low-energy 
excitations are dramatically changed. Recently, dynam- 
ics of electron spins in the field-induced critical phases 
are actively investigated by various experimental tech- 
niques^— Among these experimental techniques, elec- 
tron spin resonance (ESR) occupies a unique position 
in its sensitivity to interactions between electron spins. 
In fact, thanks to this advantage of ESR, many inter- 
esting ESR experiments have been performed in one- 
dimensional quantum spin systems under high magnetic 
field4r£ These recent precise ESR experiments enhance 
the significance of reliable quantitative theory of ESR in 
the field-induced critical phase. 

Despite the theoretical and experimental importance of 
the field-induced critical phase, ESR in the field-induced 



critical phase is less studied by theorists. This situation 
is in contrast to the fact that S = 1/2 HAF critical chain 
whose low-temperature ESR is well understood^— Note 
that ESR of the S = 1 HAF chain has been well stud- 
ied i 11 ' 12 However, these theories mainly treat the gapped 
phase of S = 1 HAF chains. Thus, it is an objective 
of the present paper to propose a theory of ESR in the 
field-induced critical phase, especially around quantum 
critical points, of one-dimensional quantum spin systems 
in an organized manner. 

In this paper, we consider an S = 1 HAF chain with a 
general form of a single-ion anisotropy 

V. = J^^Sj ■ Sj+i - g e JiBH ^2 Sj 

3 3 

+ Dj2(s*r + ey,[(s]) 2 (s-) 2 } (i) 



in the whole range of the magnetic field, from zero field 
H = to the saturation field H = H C 2- g e and (1b are 
Lande g factor of electron spin and \ib is Bohr magne- 
ton. Wc put H— ks = gefJ-B = 1 unless otherwise stated. 
Especially we focus on a shift of the resonance frequency 
(ESR shift) caused by weakly anisotropic spin-spin inter- 
actions. 

We reported, in our preceding Rapid Communica- 
tion^ that the ESR shift in a range H < H c i is well 
explained by, so called, the form factor perturbation the- 
orjsii (FFPT) around an integrable field theory. In the 
case of S — 1 HAF chain, the 0(3) nonlinear sigma model 
(NLSM) plays the role of the unperturbed integrable field 
theory in FFPT. In the Rapid Communication^ we ap- 
plied the FFPT to the analysis of the ESR shift in H 
and H as U c \ , where we utilized a close relation of effec- 
tive field theories in two different regions, H as and 
H as Hci. This paper is also intended to take a closer 
look at this remarkable feature. 
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In the next section, we will briefly review a general 
framework of perturbative treatments for the ESR shift. 
We consider ESR shifts in three regions; the low-field 
gapped region (Sec. HV| ; a region near the lower critical 
field (Sec. |VJ); and a region near the upper critical field 
(Sec. IVIj) . In each region, we introduce an effective field 
theory and apply it to analysis of the ESR shift at low 
temperature. Sec. IVIII is devoted to a comparison of our 
theory with recent ESR experiments^ of an S — 1 HAF 
compound. In App.fAJ we discuss a qualitative difference 
of the single-ion anisotropy and an exchange anisotropy 
from the viewpoint of ESR shifts. 



II. FRAMEWORK 

Here we briefly review the perturbation theory of the 
ESR shift. The ESR experiments measure an absorption 
of an oscillating magnetic field by electron spins, where 
a microwave is typically applied. From the absorption 
spectrum, we are able to extract information on dynam- 
ics of electron spins. Within the linear response theory, 
the ESR spectrum I(uj) oc wx"_(w,(|) is written in the 
retarded Green's function, 



dte lult ([S + (t),S-(0)}) 



(2) 



Here = S x ± iS v denote transverse components of the 
total spin S = £\ Sj. The correlation ([£+(*), S - (0)]} 
is SU(2) invariant. Thus, if the whole Hamiltonian pre- 
serves the SU(2) symmetry in the spin space, fl5) is triv- 
ially constant. In the presence of the magnetic field, the 
symmetry of the Hamiltonian is lowered to U(l) at most. 
If spin-spin interactions preserve the SU(2) symmetry, ([2]) 
is simple despite the presence of interactions. 



X 'l_(co,q = 0) = 2n{S*)5(oj-H). 



(3) 



The resonance frequency w r equals to the paramagnetic 
one ui r = H at any temperature. 

If spin-spin interactions do not preserve the SU(2) sym- 
metry, the above discussion breaks down and the res- 
onance frequency is shifted from the paramagnetic one. 
Let us assume that the Hamiltonian is composed of three 
terms, 



H=Hq+ Hz + H', 



(4) 



where Ho represents SU(2) symmetric interactions, Hz 
is the Zeeman term, and H' represents anisotropic inter- 
actions. The model ([T]) has the same form as ((3]). If the 
anisotropic interaction is weak, we are able to consider a 
perturbative expansion of the resonance frequency. 

The first order perturbative expansion of the resonance 
frequency was proposed first by Kanamori and Tachikii^ 
and later applied to quantum spin systems by Nagata 
and TazukeJ^±&±£ Ref. \M derived the ESR shift Suj = 



io r — H from equal-time correlations at the lowest order 
in a general formalism, 



[[H',S+],S-]) 



(5) 



The average (• ■ ■ )o is taken with respect to the unper- 
turbed Hamiltonian H^°\ 



H (0} =H + H 2 



(0) 



While we thus far treated the ESR spectrum as a func- 
tion of the frequency u) with a fixed H in above discus- 
sions, this is often not the case in actual ESR experi- 
ments. The ESR spectrum is usually obtained as a func- 
tion of H with a fixed lj. In this case, the ESR shift is 
defined as 



SH = H r - uj/g e /j, B . 



(7) 



H r is the resonance field. Note that the g factor used 
in ([7]) is determined at the high temperature limit. By 
definition, (J7J approaches zero as T — > +oo. At a low 
temperature T %J, it generaly holds that SH ^ 0. Ac- 
cording to Rcfs. l8| andl9|, within the first order accuracy, 
the ESR shift Q satisfies 



([[H',S+],S-])o 
2<S*) 



(8) 



We should emphasize that Eq. ([5]) is equivalent to ([5]). 
Therefore, as long as we are concerned with the first or- 
der perturbation theory around ((6]), it does not matter 
whether we change uj or H . 

We apply the formula ([5]) to our model (TTJ), namely, 

^^^^■■Vi-E^-^ ( 9 ) 

3 3 

H' = D J2(S^) 2 +EJ2 [(S]f (SJ) 2 } • (10) 

3 3 

The ESR shift (|5|) is, in this case, given by 

5u = f D {z)Y D (T,H), (11) 
f D (z) = D(l - 3z p 2 ) - 3E(z q 2 - z r 2 ), (12) 

YD ^ H ^^y o T,[m) 2 )o-^} (13) 

The unit vector z = H/H is parallel to the magnetic 
field. And z — (z pi z q ,z r ) is represented in the (p, g,r) 
coordinate in (p}. Thus, if we fix the orientation of the 
magnetic field, then we may regard Jd{z) as a constant. 
For simplicity, we hereafter set H = Hz where z is the 
unit vector along z axis. We call Yjj(T, H) as a normal- 
ized ESR shift. The normalized ESR shift is useful for 
our purpose because it can be applied to systems with 
any value of D and E. 
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Figure 1. (Color online) Quantum Monte Carlo results of 
the normalized ESR shift (| 13p induced by the single-ion 
anisotropy (|10[1 for temperatures T/J = 0.1 - 0.5. The sys- 
tem size is L — 40 sites. The lower critical field H c % = 0.41J 
and the upper critical field H C 2 = 4J are guided by the dot- 
ted and the dashed lines respectively. There is an extremum 
around H — H c i + T. This non-monotonic behaviour of the 
ESR shift is understood by the finite temperature crossover. 



III. QMC RESULTS 



We numerically derived the normalized ESR shift (JT3J) 
from quantum Monte Carlo (QMC) calculations. The 
QMC results of the H dependence of Yd (T, H) is shown 
in Fig. [T] We found several characteristics from Fig. [1] 
(i) The normalized shift is approximately proportional 
to H in the ranges < H < H c \ and H c \ < H < H C 2- 
The slope 8hYd(T, H) is negative in the former and pos- 
itive in the latter range, (ii) The normalized shift has a 
minimum around H = H&. The field which gives the 
minimum increases as the temperature increases, (hi) 
The normalized shift becomes zero at a certain value of 
H because Y D (T,H = 0) = and d H Y D (T, H)\ H=0 < 
hold at H = and the saturating value Yd(T, H > H^) 
is positive. Note that the field dependence in Fig. [1] is 
qualitatively different from that of S = 1/2 HAF two-leg 
ladder systems In S = 1/2 HAF two-leg ladder sys- 
tems, we fail to find the change of the sign of the ESR 
shift. Namely, the normalized shift is non-zero from the 
infinitely weak field to the saturation field (see Fig.l in 
Ref. [ljj|). The above three features suggest that the mag- 
netic field dependence of the normalized ESR shift re- 
flects the finite-temperature crossover. Below we analyze 
the normalized shift p3[) in the gapped, lower critical, 
upper critical regions. 



IV. LOW-FIELD GAPPED PHASE 
A. Effective field theory 

First we review the zero-field case, then we will extend 
the argument to the low-field case. The unperturbed 
model iJHJ) in the absence of the magnetic field, 

H = jJ2S j -S j+1 . (14) 

3 

has an excitation gap Ao = 0.41 J j2S which is called as the 
Haldane gap. Haldane proposed that HAF chains with 
an integer quantum spin number S have an excitation 
gap Ao based on a semiclassical field theory, the 0(3) 
nonlinear sigma model (NLSM)^.^ The 0(3) NLSM 
has a Lagrangian, 

C = —dun ■ d^n + —n ■ dtn x d x n. (15) 
2g 47T 

The contraction d^n ■ d^n = (d t n) 2 — (d x n) 2 was taken. 
For simplicity we put the spin-wave velocity to unity. 
The field n(t,x) represents an antiferromagnetic order: 

S x ~ y/S(S + l)(-l) x n(x) + L(x) (16) 

The uniform component L = n x dt/n/g is quadratic in 
n. The coupling constant O = 2ttS is equal to or it 
mod 27T. 

The 0(3) NLSM is integrable when 6 = 0,tt (mod 
2tt). In the case 9 = tt, the 0(3) NLSM is critical^ 
On the other hand, in the case O = of our interest, 
the 0(3) NLSM has massive triplet particles, which is 
called magnons, as the lowest excitations. The triplet 
magnons are created by n a (t,x) oc (— l) x S a (t, x) (a = 
x,y,z). Thus, the field S a satisfies the relation 

(-l) x (0|S°(t,x)|fli,oi) = S aai VZe ix ^. (17) 

|0) is the ground state, |0i,ai) is a one-magnon state 
with the rapidity 6\ and the index a\ = x, y, z, and Z is 
the renormalization factor which will be discussed later. 
The 0(3) NLSM is Lorentz invariant, and the triplet ex- 
citations obey a dispersion relation po = \J Ao 2 + pi 2 
parametrized by a single parameter 9. 

pa = Aocosh6>, p\ = Aosinh# (18) 

This parameter 9 is called as a rapidity, which uniquely 
determines the energy po and the momentum p\ of 
magnons. Therefore, the one-magnon state \9\,ai) is 
fully characterized by the rapidity 9\ and the index oi. 
We normalize the state \9,a) by 

(6 1 ,a 1 \02,a 2 )= 4*6^0,6(6! -6*). (19) 

Instead of the one-magnon state, we may consider cre- 
ation and annihilation operators of magnons, which we 
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denote Z a (6) and Z' a {9) respectively. These Z a {9) and 
Z\{ff) are called as Faddeev-Zamolodchikov (FZ) opera- 
tors and satisfy the following algebra. 

z ai (0i)z a2 (8 2 ) = s^e, - e 2 )z b2 (9 2 )z bl (e 1 ) (20) 
zKe^zKe,) = - a )i£(0 a )<(0i) (21) 

Z ai (e 1 )Zl 2 {6 2 ) = AnSa.aM - 9 2 ) 

+ s b a i a b i(e 1 -e 2 )zl(e 2 )z bl (e 1 ) 

(22) 

The factor S^f(9) is an S matrix. The S matrix possesses 
information of two-magnon scatterings. If the magnon 
cr eated by Z\(9) were a free boson (a free fermion) , the 5* 
matrix is simply S c a f(9) = 5 ad 5 bc {S c a f(9) = -6 ad S bc ). In 
reality, the magnon is neither boson nor fermion. Thus, 
the S matrix is complicated and it depends on the rapid- 
ity. Fortunately S matrix of the 0(3) NLSM is exactly 
known. 

S C J(9) = SabSdaiiO) + 5 ac 8 bd a 2 (9) + S ad S bc a 3 (9) (23) 
<7i(0)'s (i = 1, 2, 3) are 



<T 2 (9) = 



2m9 


(9 + m)(9- 


2m) 


9(9 - m 


) 


(9 + m){9- 


2m) 


2m(m' — 


0) 


(9 + iri)(9- 


2m) 



Using the FZ operator, we can create a one-magnon state 
\6,a) and its conjugate, 

\9,a) = Zt(9)\0), (6,a\ = {0\Z a (9). (24) 

n-magnon states are specified by a set of rapidities 
• • ■ ,9 n } and indices {ai, ■ ■ • , a n }j that is, 

\9 u a l ;---;9 n ,a n ) = Zl l {9 l )---Zl n {9 n )\Q), 
(0i,ai;-.. ;9 n) a n \ = (Q\Z an (9 n ) ■ ■ ■ Z ai (9x). 

They are normalized as follows: 

(9i,a\;--- ]0 n ,a n \9' 1 ,a' 1 ] ■ ■ ■ ;9 n ,a m ) 

n 

= S nm (4m) n J[8 aia[ 5{9 l -9 , l ) (25) 



1=1 



A matrix element 



F o (0i,ai;--- ;9 n ,a„) = (O|O(O)|0i, oi; • • • ;9 n> a n ) 

(26) 

is called as an n-magnon form factor of a local opera- 
tor 0(t,x). Here 0(0) is an abbreviation of 0(0,0). A 
Lorentz boost of the 0(3) NLSM alters © to 

(0|O(t,a:)|ffi,oi;-- - ;9 n ,a n ) 

= F (9 1 ,ai;--- ;9 n ,a n y( tp °- xp i\ (27) 



where Pq and Pi denote the total energy and momentum. 



P a = ^2 A cosh(9 m , Pi = ^2 A sinhfl 



m— 1 



m— 1 



(28) 



For example, the relation (jTT]) is equivalent to the one- 
magnon form factor of S a at the origin 



, ai) = (-l) r vZS a , ai 



(29) 



The relation (fl~7| connects the low-energy effective field 
theory and the physical operator S a (t,x) in the original 
spin model. The renormalization factor \/~Z inevitably 
depends on short-distance, non-universal physics and 
cannot be determined within the effective field theory. Z 
is determined only by numerical calculations. Z w 1.26 
is concluded from Density matrix renormalization group 
calculations ! 24 ' 25 It is emphasized that eq. (|29|) should 
not be interpreted as an identity between the physical 
spin operator S a and the FZ operator. The spin opera- 
tor «S a also has nonvanishing higher-order form factors. 
Thus the form factor of the powers of S a is not solely 
determined by the one magnon form factor (|29p . even in 
the leading order. 

Let us consider the tracelcss symmetric tensor 



S a6 = gagb _ ^ (30) 

T, ab has a two-magnon form factor, 
F^ab(0 1 ,ai;92,a 2 ) = -iZ 2 S ab S aia2 (3S aai - \)tp 2 (0i 



-9 2 ). 
(31) 

In the case of O(iV) NLSM^ ip 2 (6) is given by an inte- 
gral. 



MO) 

K N (x) 



sinh I — | 



00 cfa; . cosh[(7r + i9)x] — 1 
— K N (x)- 



sinh(7ra;) 



-2irx/(N-2) 



1 



(32) 
(33) 



We performed the integral and derived an explicit form 
of ^J 2 (9) for the N = 3 case in our preceding paper 



^ 2 (0) = | (0-™) taxing 



(34) 



The two-magnon form factor pip is now determined ex- 
cept for the non- universal factor Z 2 . We emphasize that 
Z 2 is an independent parameter from Z. We have de- 
termined Z 2 ~ 0.24 by using the infinite time evolving 
block decimation method^ 

The basis {\9i, a- ■ ■ ■ ; 9 n , a n )} with n = 0, 1, 2, • • • is 
complete and orthonormal. The identity 1 reads 



1 = 



OO _^ 

^1 z\2 



■d9 n 



n=l ai---a n 

x \9\,ai\ ■ ■ ■ ;0n,a n ){9i,a\; 



(47T) 



i,a n \ 



(35) 
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We note that an important relation of form factors, 
the crossing relation. In subsequent sections, we will en- 
counter matrix elements such as (9 2 , a2|0(O)|#i, ai). The 
crossing relation allows one to relate this matrix element 
to form factors. 



(0a,aa|0(O)|0i,oO 



(O|0(O)|0i,ai;0 2 -7ri J di) 
Fo(0i,ax;9 2 — 7ri,a 2 ) (36) 



The index a represents an index of an anti-magnon con- 
jugate to the magnon with the index a. If we employ 
the labeling a — x,y, z, then a = a holds. If, on the 
other hand, we employ a labeling a = +,0,—, namely 
(n+,n°,n-) =_ ((n x + i n y)/V2,n z ,(n x -in v )/V2), we 
have + = —,0 = and — = +. The latter labeling is 
required in the presence of a weak magnetic field. 

Under a weak magnetic field H < Aq , the unperturbed 
system ([9]) still have a finite gap Ao — H. Here we have 
to replace the dispersion relation (|T5| to 



(37) 



Po = Ao cosh 9 — aH, p\ = Aq sinh ( 



where a = 0, +, — . Thus, the triplet degeneracy is lifted 
by the Zeeman splitting term. If the magnetic field is 
very weak H <C Ao, then we may use the form factors 
evaluated for the H = case at the lowest order of H. 



B. ESR shift 

In the limit H,T — > 0, the density of magnons becomes 
low. It should be reasonable in this dilute limit that 
we ignore multi-magnon contributions to, for instance, 
the magnetization and the normalized shift (|13[) . We 
multiply a projection operator 



Pi 



E 



d9 

47T 



IMXM 



to an operator O so that the multi-magnon contributions 
to the average (O) are projected out. Let us consider 
O = S oo (0,a;). Using the crossing relation (|3"6"|) and the 
two-magnon form factor (|31|) . we obtain 



-iZo 



d9d9' 



M9' -9 + m)e^ Pl ^- p ^ 



/-co (4^) 2 

x (2\9,0)(9,0\-\9,+)(9',+\-\9,-)(9',-\). (38) 

Thus, in the dilute limit, the numerator J2j[3((Sj) 2 )o — 2] 
of the normalized shift (| 13[) is approximated as follows. 

£[3<(S/) 2 >o-2] 



3 / ^(S oo (0,x)) o 
-6Z 2 



-oo^E(9) \2T 



(39) 



Here E{9) = AocoshS is the zero-field dispersion. Simi- 
larly, the magnetization is given by 



(S 2 



2 sinh 



d9 

47T 



-E{8)/T 



(40) 



From (|39p and (|4TI)) . the normalized shift in the dilute 
limit reads 



Wff) = -^tanh(|; 



vd9 
AnE{9) { 



-E{6)IT 



-E(8)/T 

47T 



(41) 

PT|) correctly reproduces the features of the normalized 
shift, Yd (T, H) cx H and d H Y D (T,H) < in the limit 
H — >• 0. However, (|4"Tj) cannot explain the upturn of the 
normalized shift around H = H c \. In order to extend 
(|4"TT) to the region H ~ H c \, we must take into account 
the repulsive interaction between magnons. 



V. NEAR LOWER CRITICAL FIELD 
A. Effective field theory 

At if = A , the lowest magnon band specified by 
the index a = + touches the ground state. The point 
H c i = A corresponds to a quantum critical point. 
Above H c i, gapless excitations exist. Thus, H = H c \ 
separates the low-field gapped phase (called as the Hal- 
dane phase) and the high-field gapless phase (the field in- 
duced critical phase). We call H c \ the lower critical field. 
The quantum phase transition occurs only at T = 0. At 
finite temperature, there is a region which reflects fea- 
tures of the quantum critical point. This region is usually 
called as a quantum critical region, which is characterized 
by \H-H cl \<T. 

It is known that a free fermion theory describes low- 
energy behavior of S = 1 HAF chain in the quantum 
critical region»ii22~-2£i The free fermion has a dispersion 
relation, 



E(k) 



2A 



/'• 



(42) 



The chemical potential is fx = H — H c \. As the chem- 
ical potential of the free fermion increases, the number 
of the free fermion also increases. In terms of spin sys- 
tems, the number of the free fermion is identical to the 
magnetization density m + (T,H) = (S z ) /L. 



TA Q 
2tt 



deD(e)f(e-fi) 



Li 1/2 (- e " /T ) 



(43) 



L is the length of the spin chain, D(e) = e x l 2 is the 
density of states, and /(£) = (e^/ T + 1) _1 is the Fermi 
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Figure 2. (Color online) Comparisons of QMC and analytic 
results at (a) T/J = 0.1 and at (b) T/J = 0.2. Open sym- 
bols (circles and triangles) represents QMC data. The solid 
curves denote the normalized shift in the dilute limit (|41|l . 
The dashed curves correspond to (|47|) . 



distribution function. In the second line, the integral is 
performed explicitly, with the result given in terms of the 
polylogarithm function 

oo m 

L — ' m n 

171= 1 

Above the quantum critical region H > H c \ , a gapless 
excitation with a linear dispersion E{k) ~ k dominates 
the low-temperature physics of the 5=1 HAF chain. 
The excitation is identified with the Tomonaga-Luttinger 
(TL) liquid<2i2i We do not go into detail on the TL liquid 
in the field-induced critical phase. 



B. ESR shift 

Since Oi; • • • ; 9 n , a n }} is complete, we can expand 
the arbitrary operator 0(t, x) in the power of FZ opera- 



tors. For instance, J dxT, aa (0,x) is expanded as 
J dx T, aa {Q,x) 

= f JZ -£^[24Wo(0) - zl(0)Z + {6) 

- Zi (6)Z- {9)] + (higher-order terms). (45) 

The omitted higher-order terms contain, for instance, a 
quartic term Z\ (9 1 )Z a2 (6 2 )Zl 3 (6 3 )Z a4 (0 4 ). The projec- 
tion (|38p corresponds to an approximation which drops 
the higher-order terms of (|45|) out. To improve the re- 
sult (|4"1"T) , we need to accurately evaluate the higher-order 
terms of the expansion (|4"5|) . 

Here, without evaluating the higher-order terms di- 
rectly, we replace the FZ operators to the fermionic cre- 
ation and annihilation operators following a rule, 

Z a (9) ~ \l^c a (k), Zt(0) ~ /^ct(fc), (46) 

with k = Aq sinh 9. This rule correctly reproduces the 
anticommutation relations, 

{c a {k),c al {k')} =0, 

{4(fc),4,(fc')} = o, 

{c a (k),cl(k')} = 2nS aa ,6(k - k'). 

Note that the 5 matrix of the 0(3) NLSM j23} ap- 
proaches the fermionic 5 matrix, S^f(9) — > —S a dSbc, 
when 9 — > 0. Thus, the replacement (|4^|) is valid only 
around the bottom of th e band k y 0. This fermion has 
a dispersion E a (k) = y Aq 2 + k 2 — aH (a = 0,+,—). 
This dispersion relation is consistent with the dynamical 
exponent z = 2 of the free fermion theory (|42|) . 

The replacement (|4"6")) enables us to compute the nor- 
malized shift explicitly. 

Y D (T,H) = 322 r r Vdk - \2f (k)-f+(k)-f-(k)] 

(47) 

f a (k) = ( e (VA« 2 +t 2 -«iJ)/r + X )-i is thc Fermi distribu _ 
tion function and to(T, H) is the magnetization, 
r°° dk 

m(T,H) = J Jl[f + (k)-f-(k)]. (48) 

The analytic result (|4T]) is compared with thc QMC re- 
sults at T = 0.1J and 0.2 J in Fig. H Thc free fermion 
representation (|47p reproduces the minimum of the nor- 
malized ESR shift and, furthermore, agrees quantita- 
tively with the QMC data. 

VI. NEAR UPPER CRITICAL FIELD 

A. Effective field theory 

Thc field-induced critical phase ends at the upper crit- 
ical field H = H C 2 where H C 2 = 4J. Above the up- 
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per critical field, the spins are fully polarized, where 
the gap opens again and the low-energy excitation has 
a parabolic dispersion. Slightly below the upper critical 
field (H c i <C H — H c i < 0), almost all spins are po- 
larized. Here we may neglect the 5? = — 1 component 
anti-parallel to the magnetic field because it costs huge 
amounts of energy. Thus, the 5 = 1 spin is effectively 
described by an 5 = 1/2 spin. 



1 



si 



i 



V2 



(-iy 



(49) 



(crj, crj, <t?) is the Pauli matrices and uj = (crj ± i<rJ)/2. 
The unperturbed Hamiltonian (J9j) is transformed into an 
effective 5=1/2 XXZ chain. 



J 



1 



£E=J 
'(50) 

This is effectively written in terms of a free fermion, 



dk 
27 



with a quadratic dispersion, 



(51) 



(52) 



Here m = 1/2 J and fi = H C 2 — H arc the mass of the 
fermion and the chemical potential that the fermion feels. 
Thus, the effective theories around the upper critical field 
field and the lower critical field are isomorphic, while 
the mass and the chemical potential of the fermions are 
different. It should be also noted that the free fermion in 
each theory represents a different object with respect to 
the original spin system. 



B. ESR shift 




9000000Q 



O QMC (T/J = 0.1) 
■-- free fermion (T = 0) 
— free fermion (T/J = 0.1) 



3.5 4.0 4.5 

Magnetic field H/J 



5.0 



Figure 3. (Color online) The open circles denote the QMC 
data obtained for 40-site chains at T/J = 0.1. The solid and 
dashed curves are derived from the free fermion theory near 
H C 2- The former is T > data and the latter is T = data. 
The finite-temperature effect is irrelevant in the lower field 
range H < 3J. 



Fig. [3] In order to see the field dependence explicitly, 
we consider the T = case. The magnetization shows a 
singular dependence on the magnetic field at T = 0. 



}o = 1 V H C 2 — H 



0(H c2 - H) 



(55) 



The normalized shift at T = is shown by the dashed 
curve in Fig. |3J The free fermion theory (|51[) appears to 
work well in the entire region of Fig. [3] in the limit of T — > 
0. However, the numerical result in Fig. [T] shows a non- 
negligible temperature dependence for H < 3 J while the 
free fermion theory shows little temperature dependence. 
This corresponds the breakdown of the present picture 
based on spin flips from the saturated state, in the lower 
magnetic field. 



Using the mapping (|49[) . one can represent the normal- 
ized shift in the Pauli matrices. 

Y D (T,H) = (53) 

Here the average (• • • )o is taken by the Hamiltonian (|50|) 
of the effective 5 = 1/2 XXZ chain. A free fermion theory 
with the dynamical exponent z = 2 describes the low- 
energy physics near the upper critical field H C 2- Similarly 
to eq. (|43|) . the magnetization density (c|)o is given by 
the polylogarithm function as 

K)o = l + 2^Li 1/2 (- e (— >/*). (54) 

Substituting ([5"4"|) into ([53)1 . we obtain the explicit repre- 
sentation of the normalized shift. We show the normal- 
ized ESR shift computed by the free fermion theory in 



VII. NDMAP 

We apply our theory of ESR shifts to an S = 1 HAF 
chain compound Ni(CsHi4N2)2N3(PFg) (abbreviated to 
NDMAP^^). There are several 5 = 1 HAF chain 
compounds, for instance, Ni(C2H§N2)2(N02)C104 (ab- 
breviated to NENR2&), Ni(C 9 H24N 4 )(N02)C104 (abbre- 
viated to NTENP2&) and N^CsH^Na^NatClCU) (ab- 
breviated to NDMAZ^i). Among these 5 = 1 HAF 
chain compounds, NDMAP is most suitable to our pur- 
pose because NENP has an effective staggered magnetic 
field hJ2j(-l) j Sj and NTENP has a bond alternation 
SJ2j(~ ' Sj+i- The staggered magnetization mixes 
the singlet ground state \g) and the triplet excited states 
l e > : (e\J2 j (- 1 ) j Sj\g) ^ 0^ This mixing changes the 
selection rule of ESR. Such an interaction is uncovered 
by our theory. Though the bond alternation docs not 
induce the mixing, when H = 0, NTENP has a different 
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10 20 30 40 50 60 



Magnetic field H (T) 

Figure 4. (Color online) Comparison of the free fermion the- 
ory (|47[) with the experimental data on NDMAP at a temper- 
ature T = 0.05J.— The solid curve is the free fermion result 
(|47|l near H c \ and the solid circles denote the experimental 
data. We used the parameters J = 30.0 K, D/J = 0.25, and 
g e — 2.10. The magnetic field is applied along c axis, which 
corresponds to z p = 1, z q = z r — 0. The dashed curve is high- 
temperature paramagnetic resonance frequency uj = g^sH. 

ground state from that of Recently, NTENP has 

been field theoretically analyzed by using a sine-Gordon 
models The compound NDMAZ has very similar crys- 
tal structure to NDMAP. In fact, our theory is applicable 
to NDMAZ. But, NDMAZ has stronger exchange inter- 
action J w 70.6 K than NDMAP. The large J makes 
the experimental investigation of the field-induced criti- 
cal phase difficult because of the large H c i- 

Parameters of NDMAP are estimated as follows^ 

J«30.0K, D/J w 0.25 (56) 

The parameter E is much smaller than D. Here we con- 
sider the field orientation perpendicular to the easy plane, 
(z p , z q , z r ) = (1, 0, 0). Thus, the normalized shift is inde- 
pendent of the anisotropy E. 

Lj r = g el i B H-2DY D {T,H) (57) 

The Lande g factor is g e = 2.1l£ We substitute the free 
fermion theory near the lower critical field (jT7l) into ([5T|) 
and compare it with experimental data by Rcf. [5| (Fig.@|. 
They show semiquantitative agreement. Our theory gives 
a concrete support to the estimation ([56]) . 

Note that there is a zero point H = Hq(T) where the 
ESR shift vanishes, 

Y D (T,H o (T))=0. (58) 

In addition to the trivial solution Ho = 0, at T = 0.1 J, 
one can find a zero point Hq ~ 3 J in Fig. \5\ We show 
several cases D/J = 0.1, 0.2 and —0.1 with the fixed J. 
One will be able to experimentally observe the zero point 
H = Hq if an S = 1 HAF chain compound with smaller 
J < 15 K is found. In general, the zero point Hq(T) 




Magnetic field H/J 



Figure 5. (Color online) QMC data of the resonance frequency 
(|57p at T = 0.1J. The g factor is set to g e /J.B — 1 for simplic- 
ity. The dashed line uj = H corresponds to the paramagnetic 
resonance frequency. Several cases D/J — 0.1, 0.2 and —0.1 
are shown. Note that there is a zero point Ho ~ 3J where 
the shift (O vanishes. 

depends on the temperature T. The non-trivial solution 
Ho(T) of l|58p exists in a wide range of the temperature 
because Yd (T, H) is negative in H -c H ci and positive in 
H ~ H C 2- In contrast, as we will discuss in the Appendix, 
for the exchange anisotropy 

H' = E E J'aS?S? +1 , (59) 

j a=p,q,r 

we find that the ESR shift in the first order of the 
anisotropy does not change its sign in the entire range 
of H. 

By measuring the zero-field excitation gaps, the sym- 
metry of the Hamiltonian Q can be identified exper- 
imentally Let us suppose that the Hamiltonian has 
a uniaxial, U(l) symmetry, broken from the rotational 
SU(2) symmetry. This is consistent with a presence 
of either the single-ion anisotropy (JTUJ) or the exchange 
anisotropy ([59")) . It is usually difficult to distinguish these 
two kinds of anisotropic interactions because they often 
lead to qualitatively the same consequences in observ- 
ables. However, the presence or absence of the zero point 
Hq(T) of the shift at T ~ if cl is a clear signature which 
distinguishes the two cases. This may provide a new 
application of ESR, which posesses a high sensitivity to 
anisotropy unavailable in other types of measurements. 



VIII. SUMMARY 

We theoretically investigated the ESR shift caused by 
a weak single-ion anisotropy in the S — 1 HAF chain. We 
applied the Kanamori-Tachiki theory ([5]) to this system, 
and analyzed it both analytically and numerically. The 
formula (O is factorized to Slj = fu(z)Yn(T,H), which 
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is composed of the T, i7-indcpendcnt geometrical factor 
fn(z) and the T, if-dependent factor Y D (T 7 H). In this 
paper we call Yd (T, H ) as the normalized ESR shift be- 
cause the factor Jd{z) can be regarded as a constant if 
we fix the field orientation z. In contrast, the normalized 
shift Yd {T, H) does not depend on the field orientation. 
Thus, this factorization allows the general analysis of the 
ESR shift without specifying the parameters D and E. 

Quantum Monte Carlo calculations revealed non- 
monotonic magnetic field dependence of the normalized 
shift Yd (T, H) . The field dependence reflects the finite- 
temperature crossover of the S = 1 HAF chain, the low- 
field gapped phase (H < H c i), the field- induced crit- 
ical phase (H c i < H < H C 2). and the fully polarized 
phase (H c2 < H). We employed several effective field 
theories to explain the field dependence of Yd (T, H) in 
each phase. We used the exact form factors to com- 
pute Yd (T, H) in the dilute limit H,T -> 0. We extend 
the result in the low-field limit to the finite-field region 
H ~ H c i by replacing the FZ operators of the lowest ex- 
citations to the fermionic creation and annihilation op- 
erators. This replacement is reasonable in H < H c i and 
it works quite well (Fig. [2]) . Above H c i , the system is re- 
garded as the TL liquid. Though we did not go into detail 
of the ESR shift of the TL liquid in the field-induced crit- 
ical phase, it can be extracted from the analyses around 
H c i and H C 2- Near the upper critical field H C 2, the free 
fermion analysis is again effective (Fig. [3}. 

Our analysis is found to agree semiquantitatively with 
the experimental data of NDMAP in Ref. [ft Our theory 
correctly reproduces the approaching of the resonance 
frequency to the paramagnetic resonance frequency u> = 
gefisH. Furthermore, we predicted the existence of the 
special value Ho of the magnetic field where the ESR 
shift vanishes 8lo = 0. Such a sign change is absent in 
the case of an exchange anisotropy. 
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Appendix A: Exchange anisotropy 

We have so far considered the single-ion anisotropy as 
a resource of anisotropic interactions. In this appendix, 
we treat a perturbative exchange anisotropy instead of 
the single-ion anisotropy (fTO)) . The ESR shift caused by 
the exchange anisotropy (fBT))) is also factorized just like 




ch . 1 1 1 1 1 . 1 . 

1 2 3 4 5 

Magnetic field H/J 



Figure 6. (Color online) Quantum Monte Carlo results of 
the normalized ESR shift ()A3[) induced by the exchange 
anisotropy (|59[) for temperatures T/J — 0.1 - 0.5. The system 
size is L = 40 sites. The maximum around H — H c \ + T is 
also found in this case. The field dependence of Yj> (T, H) in 
a lower field region H < J at low temperatures T < 0.3J 
looks similar to that of Yd (T, H) except the overall sign, 
Yji(T,H) cx —Yd(T,H). In a relatively higher temperature 
T > 0.4J, the non-monotonic behavior of the normalized shift 
vanishes. 

CCD): 

5u = fj,(z)Yj,(T,H) (Al) 
fr(z)= 4(1 -3z a 2 ) (A2) 

a—p,q,r 

yAt,h) = -1—Y i Yl (3*«-i)<^s; +1 > 

\ j a=x,y,z 

(A3) 

We compute the normalized shift (|A3[) by QMC in the 
same manner as (fl~3"|) . Fig. [5] shows QMC results of 
the normalized shift (|A3|) at temperatures T/J = 0.1 

- 0.5. The normalized shift Yjr(T,H) behaves similarly 
to Y D (T, H) in a region where T < 0.3 J and H < J hold. 
On the other hand, in a higher field region H > J, the 
normalized shift quickly saturates to 1. 

First we consider the zero-field case. The effective field 
theory 0(3) NLSM works well at H = 0. When we move 
on to the continuum limit, we approximate the product 
SjSj +1 by the composite operator [S a (x)] 2 , 

S j S j+i ~ -C[S a {x)} 2 . (A4) 

The coefficient C is a non-universal constant. When 
keeping only the most relevant term S a {x)S a {x + ao) ~ 

— S(S + l)n a (x)n a (x + ao), we may assume C > be- 
cause the field n(x) is smoothly varying on x. Here ao 
is the lattice spacing and set to unity. The replacement 
(|A4|) immediately leads to Yj,(T,H) oc -Y D (T,H) in 
the infinitesimal field region H <C H c i. This relation is 
consistent with numerical results (Figs. [T] and [5]). 
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1.6 




0.2 0.4 0.6 0.8 1.0 1.2 

Magnetic field H/J 

Figure 7. (Color online) Comparisons of QMC and (|A5[) 
at T/J = 0.1 and 0.2. The open circles (T/J = 0.1) and 
triangles (T/J = 0.2) denote the QMC results. The solid 
(dashed) curve represents (KEfy with Z' 2 = 0.41 at T/J = 0.1 
(T/J = 0.2). 



Next we extend our discussion to the finite field region 
H ~ H cl in the exactly same manner with Sec. [V] We 
assume that the replacement (|A4|) is also valid under not 
so weak magnetic field H ~ H c \. Then, the normalized 
shift Yji (T, H) near H = H c i is given by 

(A5) 

We determine the phenomenological parameter Z' 2 by fit- 
ting (|A5|) with QMC data at T/J = 0.1. The fitting leads 
to Z 2 ~ 0.41. Fig. [7] shows the comparison of (|A5[) and 
QMC data for T/J = 0.1 and 0.2. The formula (gOp) 
reproduces the QMC data well. But, their agreement 
rapidly becomes worse as the temperature rises. This dis- 
crepancy stems from the saturation value Yji (T, H ) — > 1 
in the limit H —> +oo. While Yd(T, H) is negative in the 
low- field region H <C H cl , Yj>(T,H) is positive there. 
Thus, the sign change of the normalized shift does not 
occur for Yji(T,H). As we have discussed in Sec. IVII1 
this is in contrast to the behavior of Yd(T,H), which 
universally shows a sign change. 



I. Affleck, [Ph^s. Rev. B 43, 3215 (1991)] 

R. Chitra and T. Giamarchi, 

Phys. Rev. B 55, 5816 (1997)| 

M. Klanjsek, IT Mayaffre, C. Berthier, M. Hor- 
vatic, B. Chiari, O. Piovesana, P. Bouillot, C. Kol- 
lath, E. Orignac , R. Citro, and T. Giamarchi, 
|Phys. Rev. Lett. 101, 137207 (2008)] 

E. Cizmar, M. Ozerov, J. Wosnitza, B. Thielemann, 
K. W. Kramer, C. Riiegg, O. Piovesana, M. Klanjsek, 
M. Horvatic, C. Berthier, and S. A. Zvyagin, 



Phys. Rev. B 82, 054431 (2010) 



T. 
H. 



Kashiwagi, 
Miyazaki, 



M. 
I 



Hagiwara, 
Harada, 



Kimura, Z. 
and K. 



Honda, 
Kindo, 



Phys . Rev. B 79, 024403 (2009) | 
V. N. Glazkov, A. I. Smirnov^A. 



Zheludev, and B. C. 



Sales, | Phys. Rev. B 82, 184406 (2010) 
M. 



1. 



Affleck, 



Oshikawa and 
|Phys. Rev. B 65, 134410 (2002) | 
M. Brockmann, F. Gohmann, M. Karbach, A. Kliimper, 
and A. Weifie, |Phys. Rev. Lett. 107, 017202 (20lT)| 
M. Brockmann, F. Gohmann, M. Karbach, A. Kliimper, 
and A. Weifie, |Phys. Rev. B 85, 134438 (2012) | 
Y. Maeda, K. Sakai, and M. Oshikawa, 

Phys. Rev. Lett. 95, 037602 (2005)| 



I. Affleck, |Phys. Rev. B 46, 900 2 (1992)| 

L ' 1 T. Sakai and H. Shiba, |J. Phys. Soc. Jpn. 63, 867 (1994)] 

13 S. C. Furuya T. Suzuki, S. Takayoshi, Y. M aeda, and 
M. Oshikawa, [Phys. Rev. B 84, 180410 (2011)| 

14 D. Controzzi and G. Mussardo, 
|Phys. Rev. Lett 92, 021601 (2004)] 

15 J. Kanamori and M. Tachiki, J. Phys. Soc. Jpn. 17, 1384 
(1962). 

16 K. Nagata and Y. Tazuke, 



J. Phys. Soc. Jpn. 32 337 (1972)| 

K. Nagata, |J. Phys. Soc. Jpn. 40, 1209 (1976)| 



18 y 



Maeda 



and 



M. 



Oshikawa, 



J. Phys. Soc. Jpn. 74, 283 (2005)| 
S. C. Furuya, P. Bouillot, C. Kollath, M. Oshikawa , and 
T. Giamarchi, |Phys. Rev. Lett. 108, 037204 (2012)1 
S. Todo and K. Ka to,|Phys. Rev. Lett. 87, 047203 (2001) 
F. D. M. Haldane, |Phys. Rev. Lett. 50, 1153 (1983)| 
F. Haldane, Phys. Lett. A 93, 464 (1983). 
I. Affleck and F 



|Phys. Rev. B 36, 5291 (1987)| 
24 E. S. S0rensen 



Phys. Rev. B 49, 13235 (1994) | 



D. 


M. 


Haldane, 


and 


I. 


Affleck, 


and 


I. 


Affleck, 



E. S0rensen 

|Phys. Rev. B 49, 15771 (1994) | 

J. Balog and P. Weisz, Nuc l. Phys. B 778, 259 (2007). 

E. H. Lieb and W. Liniger,|Phys. Rev. 130, 1605 (1963)| 

E. H. Lie b,|Phys. Rev. 1 30, 1616 (1963) | 

I. Affleck, |Ph7s. Rev. B 41, 6697 (1990)"[ 

Y. Maeda! C~! Hotta, and M. Oshikawa, 



Phys. Rev. Lett. 99, 057205 (2007) 



R. 



M. 



Phys. Rev. B 66, 144416 (2002) | 
Z. Honda, H. Asakawa, 



Konik and P. Fendley, 

and K. Katsumata, 



Phys. Rev. Lett. 81, 2566 (1998) 



and 



A. Zheludev, Y. Chen, C. L. Broholm, Z. Hon da 
K. Katsumata, |Phys. Rev. B 63, 104410 (2001) | 
A. Zheludev, Z. Honda, C. L. Broholm, K. Katsumata 
S. M. Shapiro, A. Kolezhuk, S. Park, and Y. Qiu 



|Phys. Rev. B 68, 134438 (2003) | 

35 J. P. Renard, M. Verdaguer, L 
Erkelens, J. Rossat-Mignod, 
|Europhys. Lett. 3, 945 (1987)| 

36 A. Escuer, R. Vicent e, 
|J. Chem. Soc, Dalton Trans 

37 Z. Honda, K. Katsumata, H 
mada, T. Ohishi, T. Manabe, 



W. A. C. 

Stirling, 

Solans, 

531 (1997)| 

~K. Katori, K. Ya- 
and M. Yamashita, 



P. Regnault, 
and W. G. 

and X. 



11 



|J. Phys.: Condens. Matter 9, L83 (1997) 



T. Suzuki and S.-i. Suga, |Phys. Rev. B 72, 014434 (2005) 



Tamaki 



and 



M. 



Oshikawa, 



Phys. Rev. B 85, 134431 (2012) | 



